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ABSTRACT 


Two  controlled  source  experiments  were  conducted  in  Nevada  in  2006  and  2007  to  study  infrasound  signal 
propagation  at  distances  less  than  300  km  from  the  source.  This  area  around  a  source  is  sometimes  referred 
to  as  the  “zone  of  silence”  because,  based  on  ray  theory,  no  infrasound  energy  is  expected  to  return  to  the 
ground.  However,  many  observations  of  returned  infrasound  energy  have  been  made  at  these  distances. 

In  2006,  three  temporary  infrasound  arrays  were  deployed  at  distances  of  76,  108  and  157  from  the  source. 
In  2007,  the  site  at  157  km  was  reoccupied,  and  data  were  also  recorded  at  288  km  from  the  source.  These 
arrays  were  supplemented  with  data  from  the  permanent  Nevada  Infrasound  Array  (NVIAR),  located  37 
km  from  the  source.  In  addition,  meteorological  data  from  balloon  launched  rawinsondes  were  collected  in 
the  path  of  the  propagation. 

Interesting  results  were  derived  from  the  travel  time  analysis.  In  2006,  the  site  at  76  km  recorded  both 
tropospheric  and  stratospheric  arrivals,  while  at  108  and  157  km  only  stratospheric  arrivals  were  recorded. 
In  2007,  the  site  at  157  km  recorded  both  tropospheric  and  stratospheric  arrivals,  while  at  288  km  both 
stratospheric  and  thermospheric  arrivals  were  recorded.  While  stratospheric  arrivals  appear  most 
frequently,  atmospheric  modeling  with  the  InfraMAP  software,  did  not  predict  these  arrivals.  Current 
modeling  efforts  focus  on  the  tropospheric  arrivals. 

Amplitude  variations  at  NVIAR  for  sources  with  the  same  yield  varied  more  than  an  order  of  magnitude  on 
two  consecutive  days.  The  site  located  157  km  from  the  source  observed  variations  by  a  factor  close  to  five. 
We  therefore  attempt  to  estimate  the  yields  of  the  explosions  using  the  predominant  frequency  content  of 
the  signals.  The  physical  basis  for  such  a  method  is  found  in  an  increased  acoustic  transit  time  of  the 
explosion  blast  radius  with  increased  yield.  Past  formulas  were  developed  for  old  nuclear  atmospheric 
explosions,  but  no  extrapolation  to  the  lower  yields  was  ever  performed.  In  the  current  work  we  use 
controlled  sources  with  explosive  yields  of  2000-4,000  lbs.  in  order  to  verify  the  applicability  of  such  a 
scale  to  lower  yields.  Preliminary  results  show  that  this  is  possible. 


894 


Form  Approved 
0MB  No.  0704-0188 


Report  Documentation  Page 


Public  reporting  burden  for  the  collection  of  information  is  estimated  to  average  1  hour  per  response,  including  the  time  for  reviewing  instructions,  searching  existing  data  sources,  gathering  and 
maintaining  the  data  needed,  and  completing  and  reviewing  the  collection  of  information.  Send  comments  regarding  this  burden  estimate  or  any  other  aspect  of  this  collection  of  information, 
including  suggestions  for  reducing  this  burden,  to  Washington  Headquarters  Services,  Directorate  for  Information  Operations  and  Reports,  1215  Jefferson  Davis  Highway,  Suite  1204,  Arlington 
VA  22202-4302.  Respondents  should  be  aware  that  notwithstanding  any  other  provision  of  law,  no  person  shall  be  subject  to  a  penalty  for  failing  to  comply  with  a  collection  of  information  if  it 
does  not  display  a  currently  valid  0MB  control  number. 


1.  REPORT  DATE 

SEP  2008 

2.  REPORT  TYPE 

3.  DATES  COVERED 

00-00-2008  to  00-00-2008 

4.  TITLE  AND  SUBTITLE 

Infrasound  in  the  ’Zone  of  Silence’ 

5a.  CONTRACT  NUMBER 

5b.  GRANT  NUMBER 

5c.  PROGRAM  ELEMENT  NUMBER 

6.  AUTHOR(S) 

5d.  PROJECT  NUMBER 

5e.  TASK  NUMBER 

5f.  WORK  UNIT  NUMBER 

7.  PERFORMING  ORGANIZATION  NAME(S)  AND  ADDRESS(ES) 

Southern  Methodist  University, 6425  Boaz  Lane, Dallas, TX, 75205 

8.  PEREORMING  ORGANIZATION 

REPORT  NUMBER 

9.  SPONSORING/MONITORING  AGENCY  NAME(S)  AND  ADDRESS(ES) 

10.  SPONSOR/MONITOR’S  ACRONYM(S) 

11.  SPONSOR/MONITOR’S  REPORT 
NUMBER(S) 

12.  DISTRIBUTION/AVAILABILITY  STATEMENT 

Approved  for  public  release;  distribution  unlimited 

13.  SUPPLEMENTARY  NOTES 

Proceedings  of  the  30th  Monitoring  Research  Review:  Ground-Based  Nuclear  Explosion  Monitoring 
Technologies,  23-25  Sep  2008,  Portsmouth,  VA  sponsored  by  the  National  Nuclear  Security  Administration 
(NNSA)  and  the  Air  Force  Research  Laboratory  (AFRL) 

14.  ABSTRACT 

see  report 


15.  SUBJECT  TERMS 


16.  SECURITY  CLASSIEICATION  OE: 

17.  LIMITATION  OE 
ABSTRACT 

18.  NUMBER 
OE  PAGES 

19a.  NAME  OE 
RESPONSIBLE  PERSON 

a.  REPORT 

unclassified 

b.  ABSTRACT 

unclassified 

c.  THIS  PAGE 

unclassified 

Same  as 
Report  (SAR) 

8 

Standard  Form  298  (Rev.  8-98) 

Prescribed  by  ANSI  Std  Z39-18 


2008  Monitoring  Research  Review:  Ground-Based  Nuclear  Explosion  Monitoring  Technologies 


OBJECTIVES 

The  primary  objective  of  this  research  study  is  the  propagation  of  infrasonic  waves  at  distances  up  to 
300  km  from  the  source.  This  area,  extending  up  to  300  km  from  the  source  has  been  called  the  “Zone  of 
Silence”  because  classical  ray  theory  does  not  predict  any  infrasound  arrivals.  However,  under  appropriate 
wind  and  temperature  conditions,  infrasound  recordings  in  the  “Zone  of  Silence”  were  reported  several 
times  (McKenna,  2005;  Golden  et.  al.  2007). 

RESEARCH  ACCOMPLISHED 

Two  controlled  source  experiments  were  carried  out  in  2006  and  2007  near  Hawthorne,  NV,  where 
ammunition  disposal  operations  take  place  on  a  regular  basis.  Shown  in  figure  1  is  the  location  of  the 
Nevada  Seismic  Array  (NVAR)  and  the  configuration  of  infrasound  array.  Also  shown  is  the  location  of  the 
disposal  pits  (dubbed  New  Bomb)  and  NVAR  recordings  of  a  typical  disposal  operation.  In  a  typical 
operation  day  a  suite  of  five  explosions  is  recorded  at  NVAR,  but  the  number  of  explosions  in  a  day  could 
vary  from  3  to  10,  depending  on  the  nature  of  operation. 


Figure  1.  Location  of  NVAR  and  the  ammunition  disposal  site  (New  Bomb).  Also  shown  are  the 
NVIAR  configuration  and  a  typical  NVIAR  recording  of  the  detonations. 

During  the  2006  experiment,  three  4-element  temporary  arrays  were  installed  north  of  New  Bomb,  at 
distances  of  76,  108  and  157  km.  The  array  at  76  km  was  located  within  the  village  of  Shurz  NV  (SHURZ), 
while  the  other  2  arrays  were  north  and  south  of  the  city  of  Fallon,  NV,  designated  Fallon  north  (FALN) 
and  Fallon  south  (FALS),  along  State  Highway  95.  In  the  2007  experiment,  the  site  at  FALN  (157  km)  was 
reoccupied,  and  a  new  temporary  array  was  installed  at  288  km  from  the  source,  north  of  Gerlach,  NV 
(GERL).  The  locations  of  the  temporary  arrays  and  the  individual  layouts  are  shown  in  Figure  2.  Each  of 
the  portable  arrays  consisted  of  4  IML  infrasound  sensors  located  in  line  with  the  direction  of  propagation 
of  the  signals  at  separation  distances  of  50,  100,  and  150  meters.  The  linear  configuration  was  chosen  to 
maximize  the  time  delays  across  the  array.  In  the  2006  experiment,  data  were  collected  for  three  days, 
while  in  the  2007  experiment,  data  was  collected  for  four  days. 

In  addition  to  the  infrasound  array  data  we  collected  meteorological  data  at  the  Hawthorne  NV  airport.  We 
employed  the  National  Severe  Storms  Laboratory  (NSSL)  to  collect  balloon  launched  rawinsonde 
meteorological  data  during  our  three  experiment  days.  The  rawinsondes  were  used  to  collect  temperature, 
dewpoint,  barometric  pressure,  GPS  latitude,  longitude,  altitude,  and  wind  speed  and  direction  data. 
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Figure  2.  Satellite  image  of  the  area  of  temporary  infrasound  arrays  deployed  for  the  experiments 
with  plan  layouts  of  each  array. 
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Figure  3.  Infrasound  mean  travel  velocity  (celerity)  at  each  array  for  all  observed  arrivals. 


Interesting  results  were  derived  from  the  analysis  of  the  travel  time.  Celerity  (distance  divided  by  the  total 
travel  time)  was  used  as  an  indicator  for  the  nature  of  recorded  signals  (Muntschlecner  and  Whiteker,  1999; 
Kulichkov  et  ah,  2000).  The  celerity  values  are  shown  in  Figure  3.  In  2006,  most  of  the  recorded  signals 
were  stratospheric  (celerities  of  310-284  meters/second).  Only  SHURZ  (76  km  form  the  source)  recorded  a 
low  tropospheric  arrival  (Figure  4)  for  the  first  two  days  of  the  2006  experiment  (celerities  around 
340  meters/second).  In  the  2007  tropospheric,  arrivals  were  recorded  at  FALN  during  all  days  of  the 
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Figure  4.  Data  recorded  at  SHURZ  (76  km)  during  the  2006  experiment.  Approximately  30  seconds 
before  the  arrival  of  the  main  wavetrain  a  low  tropospheric  arrival  was  recorded.  Any 
other  subsequent  tropospheric  arrivals  would  be  masked  by  the  much  higher  amplitude 
stratospheric  signals. 


Figure  5.  Tropospheric  and  stratospheric  signals  recorded  at  FALN  (157  km)  during  the  2007 

experiment.  After  the  tropospheric  signals  from  the  first  two  explosions,  there  is  mixture  of 
stratospheric  and  tropospheric  signals. 
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experiment.  Figure  5  shows  the  FALN  recording  for  Julian  day  255.  After  the  first  two  tropospheric 
signals,  there  is  a  wavetrain  that  is  interpreted  as  a  mixture  of  stratospheric  and  tropospheric  signals.  This 
combination  of  tropospheric  and  stratospheric  signals  was  observed  for  the  first  three  days  of  the  2007 
experiment,  while  in  the  fourth  day  only  tropospheric  signals  were  observed.  For  the  array  at  Gerlach 
(288  km),  signals  were  observed  for  only  two  days  of  the  experiment,  while  during  the  other  two  days  of 
the  experiment,  the  local  noise  conditions  were  significantly  higher.  GERL  recorded  a  combination  of 
stratospheric  and  thermospheric  arrivals  (celerities  of  296  for  stratospheric  and  240  for  thermospheric 
arrivals).  The  first  thermospheric  signal  is  recorded  before  the  arrival  of  the  last  stratospheric  signals 
(Figure  6).  No  tropospheric  arrivals  were  observed  at  GERL. 


Figure  6.  Thermospheric  and  stratospheric  arrivals  recorded  at  GERL  (288  km  from  the  source). 

The  first  thermospheric  arrival  is  observed  after  four  stratospheric  arrivals,  before  the  last 
stratospheric  signal. 

In  addition,  the  carefully  controlled  explosive  yields  enable  us  to  investigate  various  scaling  relationships 
for  yield.  As  we  have  observed  in  the  past,  amplitude  variations  more  than  one  order  of  magnitude  for 
infrasound  signals  recorded  in  similar  conditions,  we  decided  to  use  a  relationship  between  the  yield  and 
dominant  period  of  the  signal.  Such  a  relationship  was  developed  from  a  suite  of  atmospheric  nuclear 
explosions: 

7  =  2.38-r^-^^  (1) 

Spectral  estimates  are  typically  obtained  using  Fourier-based  methods  which  provide  a  relative  measure  of 
the  energy  at  each  frequency.  For  simple  periodic  signals  the  estimates  might  be  accurate,  but  for  real  data 
with  noise  the  Fourier  methods  provide  a  very  poor  estimate  of  the  spectrum.  Various  windowing 
techniques  were  developed  to  reduce  the  variance  of  the  spectral  estimates  (Blackman  and  Tuckey,  1959; 
Welch,  1967),  but  they  tend  to  spread  the  energy  over  adjacent  frequencies.  An  alternative  way  to  obtain 
spectral  estimates  is  the  AR  method  (Priestly,  1981).  The  AR  method  is  a  parametric  method,  widely  used 
in  statistics  and  has  direct  applications  in  many  areas  of  interest.  The  method  provides  an  estimate  of  the 
spectrum  and  the  fundamental  (or  system)  frequencies  of  the  time  series.  For  a  quick  comparison  of  the 
methods  Figure  13  shows  the  power  spectrum  estimate  of  a  time  series  obtained  with  a  simple 
periodogram,  a  windowing  technique  (Welch,  1967),  and  an  AR  method.  By  far  the  best  estimator  is  the 
AR  method. 

An  important  problem  in  using  AR  methods  is  the  order  selection.  The  higher  the  order,  the  better  the 
spectrum  estimates  will  be,  but  the  more  difficult  it  is  to  interpret.  In  our  case,  the  order  selection  was 
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performed  empirically  so  that  the  potential  noise  bias  on  the  frequency  estimate  is  reduced.  By  adding  noise 
before  the  arrival  of  the  infrasound  signal,  and  estimating  the  spectrum  of  signal  and  noise  repeatedly,  we 
found  out  that  for  orders  higher  than  1 6  the  noise  bias  is  significantly  removed,  even  in  the  case  of  a  signal 
with  a  low  signal  to  noise  ratio.  The  accuracy  of  the  frequency  estimates  with  an  AR(16),  for  which  the 
coefficients  were  determined  with  Burg’s  maximum  entropy  method  (Burg,  1972),  was  tested  with  a 
bootstrap  algorithm  and  it  was  found  to  be  extremely  reliable.  The  standard  deviation  of  5,000  realizations 
was  determined  to  be  below  0.15  Hz  (Figure  7). 


Frequency  (Hertz) 

Figure  7.  Comparison  of  AR  and  Fourier  methods. 
CONCLUSIONS  AND  RECOMMENDATIONS 


Infrasound  modeling  was  performed  with  the  InfraMAP  software  developed  by  BBN  Technologies.  Shown 
in  figure  8  are  the  results  of  the  modeling  using  the  Naval  Research  Lab  Ground-to- Space  (NRL  G2S) 
model  for  the  2007  data.  The  G2S  model  is  supposed  to  include  the  local  meteorological  variations.  The 
left  plot  shows  the  ray-tracing  results,  and  the  right  shows  the  pressure  amplitude  field  (in  dB)  obtained 
with  the  parabolic  equation  (PE)  code.  The  locations  of  our  arrays  are  shown  on  the  left  plot.  On  that 
particular  day  FALN  recorded  tropospheric  and  stratospheric  signals,  while  at  GERL  we  observed 
stratospheric  and  thermospheric  signals.  None  of  the  actual  arrivals  are  predicted  with  either  of  the  method. 
Thermospheric  rays  are  predicted  to  bounce  at  more  than  350  km  away  from  the  source,  while  raytracing 
shows  no  stratospheric  or  tropospheric  rays.  These  results  are  similar  to  the  PE  code. 

Current  modeling  efforts  focus  on  using  the  actual  meteorological  data  acquired  in  the  paths  of  the 
propagating  signals  to  determine  the  presence  of  any  potential  duct  in  the  troposphere  that  could  have 
allowed  the  transmission  of  energy  to  greater  distances.  In  addition,  this  controlled  dataset  will  allow  us  to 
investigate  the  relationship  between  the  celerity  and  the  phase  velocity  (which  in  ray  theory  is  supposed  to 
be  equal  to  the  effective  sound  speed  velocity  in  the  place  where  that  particular  ray  is  turning). 
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Figure  8.  Ray  tracing  (left)  and  PE  calculations  through  the  NRL  G2S  model  for  the  2007  data.  The 
model  is  supposed  to  include  the  meteorological  effects  below  20  km. 
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Figure  9.  AR  estimates  for  different  days  for  2006  and  2007 
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Frequency  estimates  for  2006  and  2007  datasets  were  obtained  with  an  AR(16)  method.  Figure  9  shows  the 
Power  Spectral  Density  Functions  for  days  in  which  the  amount  of  ordnance  disposed  were  the  same  in 
total  weight  (3,809  lbs.).  However,  the  spectral  estimates  are  different.  The  mean  frequency  estimate  for  the 
2006  dataset  (40  samples  in  total)  is  0.92  Hz,  while  the  mean  for  the  2007  (40  samples)  dataset  is  1.18.  This 
difference  is  much  larger  than  the  error  of  the  AR  method,  and  therefore  they  can  be  attributed  only  to 
differences  in  the  source.  The  logs  from  the  disposal  contractors  indicate  that  though  the  total  weight  of  the 
explosives  was  similar  in  the  two  years,  the  materials  have  different  chemical  characteristics  (explosivity). 
In  2006,  the  total  of  3,809  lbs.  was  composed  of  1,800  lbs.  HEX  and  2,009-lb.  grenades,  while  in  2007,  the 
composition  was  450  lbs.  HEX,  1350  lbs.  composition  E  and  2,009  lb. -grenades.  HEX  is  the  most 
energetic  in  this  combination,  and  we  believe  that  the  higher  frequency  estimated  for  2007  is  due  to  less 
HEX.  The  operations  at  Hawthorne  will  be  on  going,  and  in  the  future  there  are  scheduled  detonations  of 
HEX  only.  We  will  use  these  HEX  only  detonations  as  references  for  future  investigations  on  the 
yield/dominant  period  relationship. 
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